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Abstract 



o 



iWe suggest baryon-baryon correlations as an experimentally accessible signature for a first-order phase transition between 
'a baryon-rich phase, like quarkyonic, and a baryon-suppressed hadronic phase in the QCD phase diagram. We examine 
the consequences of baryon-rich bubble formation in an expanding medium and show how the two-particle correlations 
Vary in the transverse and longitudinal direction depending on the strength of the radial flow, the bubble temperature, 
and the time when the baryons are emitted. 
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1,1. Introduction and Motivation 

Exploring the Quantum-Chromodynamics (QCD) phase 
diagram has been the focus of significant experimental and 
theoretical research in the last few decades. In particu- 
lar, the effect of heating and compressing nuclear matter 
'on confinement and chiral symmetry breaking has been 
studied, new phases of matter have been suggested, and 
"the existence of a critical endpoint. predicted Rel- 
ativistic heavy-ion experiments also aim to explore the 
temperature-baryon density plane and to provide evidence 
for a critical endpoint or a first order phase transition. 
Therefore, it is important to define reliable signatures for 
identifying a critical point, a new phase, or a first-order 
transition line. 

In attemtping to understand QCD, it has often proven 
instructional to explore QCD-like theories in which some of 
the parameters of the Lagrangian (quark mass rriq, num- 
ber of colors TVc, number of flavors Nf) are varied Q. 
McLerran and Pisarski used such a study to argue that 
cit large and sufficiently small temperature to baryon 
chemical potential ratio T / the phase diagram exhibits 
a strongly first-orderphase transition to what they termed 
quarkyonic matter [2l|. For large Nc and all TV/, baryon 
number has been suggested as an order parameter for the 
transition In large the baryon mass is of the or- 
der Mb ~ NcAqcd- For moderate temperatures, T ~ 
Aqcd, the expectation value of baryon number is rts ~ 
exp{fiB/T - Mb/T) ~ e-^=. This is negligibly small 
(identically zero at large Nc), and it remains that way 



as long as /is < Mb- For larger ^b the baryon num- 
ber density becomes non-zero. In the deconfined quark- 
gluon plasma phase there are no baryon masses, so that 
there is no baryon- number suppression Figure [T] shows 
the phase diagram with a hadronic (niesonic) phase and a 
high energy density quarkyonic phase. At small tempera- 
tures the transition from the mesonic to quarkyonic phase 
should be first order Q. Although large Nc is not QCD, it 
is not excluded that the QCD phase diagram bears some 
resemblance to the large Nc diagram. For realistic Nc and 
Nf the boundary of these regions may be crossovers. 
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Figure 1: Possible phase diagram of QCD. 

The question we asked ourselves was, if there is a first- 
order transition between the quarkyonic and mesonic phase, 
what are the observable consequences? The arguments 
laid down in the rest of the paper are, however, more gen- 
eral, because they only require the existence of a baryon- 
rich and a baryon-poor phase. Phase conversion between 
these can happen either by bubble nucleation or by spin- 
odal decomposition. Here we explore the phenomenolog- 
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ical consequences of bubble nucleation in QCD matter. 
Spinodal decomposition has been addressed in It has 
been suggested by Voloshin ^ that correlations in coordinate- 
space can be transferred to momentum-space through ra- 
dial flow. Baryon number fluctuations were discussed in 

01 , but the effects of radial expansion and azimuthal corre- 
lations were not considered in that work. We focus on the 
transfer of spatial correlations between baryons confined 
to baryon-rich regions, into baryon-baryon correlations in 
momentum space via radial flow. 

2. Baryon-rich Bubbles 

In this work we consider a bubble nucleation scenario, 
similar to the one proposed for strangelet condensation in 
the early universe Imagine an expanding system in 
which the initial temperature is higher than that of the 
transition between quarkyonic and mesonic matter. Ini- 
tially there is non-zero baryon number density. As the 
system cools a small amount below the phase transition 
temperature, bubbles of stable mesonic phase start to nu- 
cleate. It has been shown in a dynamical calculation by 
Paech and Dumitru Q that baryon inhomogeneities do 
form. The mesonic bubbles expand rapidly, reheating the 
system back to the transition temperature. The system 
then is in a mixed phase, which at first is mostly quarky- 
onic with small droplets of liadron gas. Then these mesonic 
droplets increase in size. At the time when roughly half 
the matter is in the quarkyonic phase one can think of 
bubbles of quarkyonic matter evaporating in the hadron 
gas, as shown in Figure [2l The ratio of baryon densities 
of the hadronic phase and that of the quarkyonic phase 
n-n/nQy ~ e'^" is tiny (identically zero in large Nc). Be- 
cause of this exponential suppression of the baryon density 
in the mesonic phase, all of the baryon number remains 
trapped inside the bubbles of quarkyonic phase, and at 
later times it can concentrate in the evaporating droplets 
of quarkyonic matter. 




Figure 2: Evaporating quarkyonic bubbles Qy in the hadron gas H. 

If bubbles of quarkyonic matter have baryon number 
larger than one, large variation will exist in the spatial 
baryon number-density. In the following we study how 
inhomogeneous the corresponding distribution in momen- 
tum space will be. The longitudinal position inhomo- 



geneities will be translated into inhomogeneities in rapidity 
because of the familiar correlation between momentum- 
and coordinate-space rapidity. In the transverse space 
there should be a similar correlation due to transverse flow 
0. 

Baryon correlations due to transverse flow are not likely 
to be washed out by thermal noise. The momentum im- 
parted to a hadron due to flow depends on the mass and 
flow velocity pfiow Mvfiow while the thermal noise de- 
pends on the mass and temperature Pnoise ^ V MT. Since 
Pfiow/Pnoise ^ Vfiow^/M/T and since the baryon mass 
is at least a factor of five larger than the temperature, 
we expect space inhomogeneities for baryon number to 
be observable as baryon-baryon correlations in momen- 
tum space. We investigate these by constructin g a Bubble 
Monte-Carlo Blast- Wave Model similar to Ref. |10j. 

3. Bubble Monte-Carlo Blast- Wave Model 

Our model has two components: One is the Bubble 
Monte-Carlo, in which the spatial distribution of baryon 
number and bubbles is generated; The other is a Monte- 
Carlo Blast- Wave which provides the momentum distribu- 
tion of baryons. We study how spatial correlations in the 
source function are translated into momentum space and 
examine the effect of baryonic bubble formation on the 
two-particle correlation functions. 

In our Bubble Monte-Carlo, we first make an assump- 
tion for the average number of baryons per bubble and 
the distribution of the baryon number per bubble. We 
generate baryon number Nb for each event according to a 
gaussian distribution 

dN/dNB^exp{~{NB-NBf/2NB}, (1) 

in a finite window of rapidity within which we know the av- 
erage baryon number, Nb- We then distribute the baryons 
into bubbles with the number of baryons per bubble ub 
taking the form 

dN / duB ^ exp{—{nB ^ tlb)'^ I'^.tlb} ■ (2) 

We sample this baryon-per-bubble distribution to gener- 
ate bubbles until we have generated the total number of 
baryons within the given window. This fixes the number of 
bubbles. Our treatment of Nb and ris is somewhat arbi- 
trary, but in this analysis rather than focusing on number 
fluctuations sensitive to this implementation, we will study 
the shape of the two particle correlations in rapidity and 
azimuth {(f). 

The bubbles are then distributed in space with two 
different transverse density profiles: The Center Weighted 
Source profile, assumes a bubble density given by a Woods- 
Saxon distribution integrated along the longitudinal z-direction. 
This distribution yields a density largest at the center of 
the system. The Surface Weighted Source gives a prefer- 
ence to bubble formation near a surface, motivated by the 
picture of an expanding system where bubble nucleation 
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Figure 3: The source functions with emission weighted most heavily 
from the center (left) or from the surface (right). 



happens preferentially at a freezeout surface, near the edge 
of the system. Transverse projections of both source func- 
tions are shown in Figure [31 For both cases, we take the 
distribution to be flat in z. Baryons are then distributed 
inside the bubbles using again a Woods-Saxon distribution 
of the form 

dN/dxdydz - 1/(1 +exp{-(rfc - R)/a}) . (3) 

The bubble radius is taken to be — 0.6 fm(nB)^/'^, so 
that the volume per baryon remains constant, the radius 
of the system R= 12 fm, and the width of the distribution 
a — rb/8 . In this way we obtain the positions (x, y, z) for 
each baryon. This discrete set of coordinates is later taken 
as the source function $7. Figure |4] shows the baryon- and 
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Figure 4: Transverse spatial baryon- and bubble-distribution for a 
center weighted source. The circles representing the bubbles are 
drawn with a fixed radius that does not reflect the variable bubble 
size r;,. 

bubble-distribution in the transverse plane of one event. 

Each bubble is then assigned a radial boost velocity as 
is done in the blast-wave model. To construct a blast-wave 
Monte-Carlo we calculate the emission function, S. The 
emission function is the probability of emitting a baryon 
of transverse mass = -y^p^ -f m-^, transverse momenta 
Pt, and rapidity Ye, at an azimuthal angle (pB from a 
boosted source of VI at temperature T. We use the discrete 
set of points generated with our bubble Monte Carlo as our 



source function VL. Following flTl, we write the emission 
function as 



S — wlt cosh(77 — YB)^{r, (ps)' 



,-(r-ro)V2Ar" 



where: 



ri=l 



— sinh/3(r, and 6 = — cosh/3(r, (^i^) . (5) 



Here we will consider central collisions, which implies the 
direction of the boost (3 is identical to the spatial azimuthal 
angle <j)b — (j^s- Generalization to finite impact param- 
eter presents no major technical difficulties. We allow 
the bubbles to decay isotropically in their flow boosted 
rest frame. Assumption of boost invariance manifests in 
the equality of the longitudinal flow rapidity to the space- 
time rapidity rj. For simplicity, we assume that the emis- 
sion is instantaneous and replace the Gaussian in ^ with 
^{t — tq). We are left to specify the transverse rapidity 
boost (3 = r[/?o + /32 cos(20f,)]. The parameter /3o is the 
transverse rapidity in the outward (ps direction. We take 
(52 to be zero. We will investigate different values for (5^ 
taking guidance from the flts in illj, which reproduces well 
the experimental data. 

4. Baryon-Baryon Correlation Results 

We investigate how baryon-baryon correlations in terms 
of relative angle (azimuth A0 and rapidity Ay) vary with 
the freeze-out temperature T, average flow velocity (/3), 
number of baryons per event Nb^ number of baryons per 
bubble ub, time at which bubbles decay Tdecay, and source 
function (Surface vs Center Weighted). We present the 
correlations in terms of a per particle correlation measure 
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Ap _ d^N/{d(j)id(j)2) - {dN/d(j)i) ■ {dN/d(j)2 



^{dN/d<P^) ■ {dN/dc^2) 



(6) 



Here p is the pair density and pref the product of single 
particle densities. This normalization of the correlation 
measure means that for all other parameters kept flxed 
Ap/^Pref is independent of the total number of baryons 



per event Nb 12[. This is because normalizing by ^ Prej 
instead of prej removes the dilution factor of one over the 
number of particles analyzed. This was confirmed to be 
true in our Monte-Carlo (the variation of rf, with Nb does 
not seem to infiuence the correlation function) . We gener- 
ate our reference spectrum by sampling the single particle 
dN/dij) and dN/dy distributions generated in our Monte- 
Carlo enough times so that the same number of pairs are 
created in our reference distribution as in our real distribu- 
tion. This ensures that Ap/ ^p^ef will integrate to zero. 
For this reason we are insensitive to a global offset and 
focus instead on the shape of the correlation, which is eas- 
ier to measure experimentally. The variable Ap/ ^ p^ef as 
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Figure 5: Proton-proton correlations in terms of azimuthal angle between them for different model parameters. For details see the text. 



a function of the azimuthal angle between baryons from 
the decayed bubble for a variety of different model inputs 
is shown in Figure [5] Bubble nucleation and radial flow 
leads to small angle correlations. We find that all corre- 
lation shapes can be reasonably well parameterized by a 
Gaussian centered at Acj) = and a constant offset. The 
correlation can therefore be summarized in terms of the 
Gaussian amplitude A and root mean squared width cr. 
By construction the quantity Ap/^/Pref is independent of 
the total number of baryons per event 12| . We confirmed 
this in our simulation, and then fixed Ng = 100. 

The top left panel in Figure [5] shows results with a 
Surface Source and T = 100 MeV, (/3) = 0.75, and ng = 
2, 4, 8. The width of the correlation is narrowest for ns — 
2, since this corresponds to the smallest bubble size, but 
the variation of a with ub is very weak (only changing from 
0.487 for nij = 2 to 0.499 for = 8). The dominant effect 
of varying hb is on the amplitude of the correlation, which 
changes from 1.3 to 8.9, respectively. This dependence 
can be understood in the following way: The number of 
pairs that can be formed from one bubble is nB{nB — 
l)/2. So the number of correlated pairs in an event is 
nsinB — l)/2 times the number of bubbles, which is the 
number of baryons Nb divided by the number of baryons 
per bubble ub- Thus the number of correlated pairs in an 
event is Nb (ns — 1)/2. Recall, that the variable Ap/ ^/p^ 
measures the number of correlated pairs per particle. For 
our simulation cx ^^("^^--1) — ^„ _ Thus the 



/Pre/ 



correlation amplitude scales with — 1 ■ Accordingly, the 



riB — 8 amplitude is approximately 7 times larger than for 
hb = 2. We note here that this amplitude only applies to 
the case that the correlation is formed only with baryons 
produced from bubbles. Any background contribution (for 
instance from baryons produced outside of bubbles) will 
dilute the signal. 

In the top right panel in Figure [5] we show the depen- 
dence of the correlation on the temperature of the bubbles 
emitting the baryons, for a fixed (13) = 0.75 and ub = 8. 
This temperature is not necessarily the temperature at 
which the bulk of the hadrons freeze-out, but could be 
closer to the value of the critical temperature Tc, where 
the bubbles actually form. We vary this T from 120 MeV 
to 180 MeV. It is clear that the shape of the correlation 
changes appreciably with temperature: higher T lead to 
a broader correlation function. In this temperature range 
and for this fiow velocity the r.m.s. width is given by 
the fit cr = 0.694 - 0.00639T + 0.0000435r2. Since the 
random thermal boost competes with the fiow boost, a 
larger bubble decay temperature washes out the correla- 
tion. Therefore, if the first order phase transition is at a 
high Tc, correlations from bubbles will be harder to detect. 

Since we expect a competing effect of the radial flow 
and the temperature, we now look at how flow influences 
the correlation. We illustrate this for T — 100 MeV, 
UB = 8, and (/3) = 0.75,0.55,0.40 on the bottom left 
panel of Figure [5] The correlation is strongly dependent 
of the mean flow velocity and is the strongest for bub- 
bles with larger boost. Remember, that flow is required to 
convert coordinate-space correlations to momentum space. 
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So again, if the phase transition happens early, before flow 
builds up, bubbles will be hard to detect. 

The bottom right panel of Figure [S] displays the flow 
effects at T = 100 and ti^ = 8, with either a Center- 
or Surface Weighted source function. The mean radial 
position of the source is larger for the Surface Weighted 
source: 10.1 fm compared to 7.4 fm. This leads to a 
larger average boost velocity: (/3) — 0.75 compared to 
= 0.55 , respectively. A larger leads to a more 
pronounced correlation structure with the Gaussian nar- 
rowing from (7 — 0.61 to ct = 0.50. We also compared a 
Surface Weighted source with a reduced /3o such that the 
(/3) is matched to the Center Weighted case. The Surface 
Weighted correlation with (/3) — 0.55 is only slightly nar- 
rower than the Center Weighted correlation with the same 
(/3). The primary difference between Surface- and Center 
Weighting is caused by the change in {(3), but some devi- 
ation occurs due to the difference in geometry. Our inter- 
pretation is that the Surface Weighted source has a smaller 
range of /? values, while the Center Weighted source has a 
broad distribution of /3 values. In this case, even for the 
same {(3) value, the Center Weighted source gives a convo- 
lution of a broad range of /3 values which leads to a broader 
correlation function. Also, a bubble near the center cov- 
ers a larger solid angle. We investigated three {(3) values 
for each geometry. The widths within the ranges studied 
can be parameterized as follows: a = 1.317— 1.103/3 for 
the Center Weighted source (see lower right panel) and 
a = 1.357 — 1.367/3 for the Surface Weighted source (see 
lower left panel). These parameterizations of our model 
are vahd for 0.29 < /3 < 0.75. 




Figure 6: Proton-proton correlation versus azimuthal angle and ra- 
pidity difference between the protons for bubble decay at times times 
Tdecay = 1 fm (left) and Tdecay = 9 fm (right). 

We also looked at correlations between baryons in the 
rapidity y direction. FigurejBlshows the proton-proton cor- 
relation as a function of rapidity difference and azimuthal 
angle between the protons at different proper times of bub- 



ble decay, i.e. Tde 



1 fm (left panel) and Tde 



9 fm 



(right panel). We found that if the time at which bubbles 
decay Tdecay is long then narrow correlations can build 
up both in A(/) and Ay directions. These correlations are 
increased in amplitude and narrowed in width. We con- 
clude that if bubble decay happens near the beginning of 



the system evolution (small Tdecay) implying a high decay 
temperature and low flow, baryon-baryon correlations will 
be broad in both directions. For bubble decay near freeze- 
out (large Tdecay), ou the other hand, strong correlations 
with narrow widths in both directions can be expected. 
Our fit for the change in the width in rapidity direction 
with time is ay = 0.37 -I- QA7t~'^-^^. 

Our studies of Ap/y^p^e/ show that, the amplitude and 
width of the correlation varies substantially within this 
bubble Monte-Carlo blast-wave model depending on the 
parameters assumed: We have shown in Figure [S] That 
correlation caused by baryonic bubbles emitted from a 
boosted source can change by a factor of two even with rea- 
sonable values for the average boost velocity and freeze-out 
temperature. The signature of a first order phase transi- 
tion will be strong if the transition happens near the freeze 
out since then the flow will be largest and the temperature 
lowest. We note, that the decay time is linked to the tem- 
perature and flow, and with a dynamical model one could 
make this connection explicit. 

We also consider how this correlation can affect mea- 
surements of elliptic flow V2 [l^ and it's fluctuations [l^ . 
One method for estimating V2 = (cos(2(/)— ^E*)) when the re- 
action plane angle 5* is not known is to calculate (cos(2A(/i)) 
from two-particle correlations. This quantity depends on 
(u2)^, cr^a' ^^'^ other correlations not related to the reac- 
tion plane, called non-flow 82- Since we only consider cen- 
tral collisions in this study (W2) is zero. The calculation 
of (cos(2A(/))) is therefore related to either non-flow corre- 
lations or V2 fluctuations. In either case, the correlations 
from our model will lead to a deviation between the 2"^^ 
and 4*'* order cumulants [l3| such that U2{2}-^ — U2{4}^ is 
non-zero. Since V2 is zero the expected deviation between 
the cumulants is given by 



V2{2Y'V2{^Y 



dA<j)cos{2A(l))dnpairs/dA(f> 
JJ" dA(l)dnpairs /dAcf) 



(7) 



This quantity depends strongly on the variables explored 
in the discussion above, as well as on the number of un- 
correlated particles included in the analysis (misidentifed 
mesons for example). Larger multiplicities will dilute the 
correlation leading to smaller differences between V2{2}^ 
and W2{4}^. The quantity A^b(w2{2}2 - W2{4}^), however, 
should be independent of Nb if all other paremeters are 
fixed because. The factor of Nb cancels out the dilution 
of W2{2}^ — W2{4}^ arising from the increase of the num- 
ber of possible pairs when the multiplicity increases. To 
confirm this we kept other parameters fixed and varied 
Nb from 60 to 200. The resulting Nb * (w2{2}^ - W2{4}^) 
did not change. With T = 100 MeV and (^) = 0.55, we 
found that Nb{v2{2}^-V2{4.}^) « 5 for = 8 and « 2.7 
for ns = 4. Increasing the temperature to 180 MeV with 
= 4 reduces Wb(-U2{2}^- W2{4}^) to « 0.37. The intro- 
duction of a number M of uncorrelated background par- 
ticles (pious for example), reduces Nb{v2{2}'^ — f2{4}^) 
by a factor of M^. Putting everything together we find 
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W2{2}^ — W2{4}^ = NbM^ ' where C can be expected to be 
between 0.3 and 5. 

Measurements of p-p and p-A correlations have been 
carried out at energies ranging from ^/s ^ 2 GeV to 200 
GeV 15| . In these HBT analyses correlation functions are 
defined as C{k*) = N pair, reai/N pair, mixed- The variable 
k* = Qinv/'^ is the relative momentum of the particles 
in the pair rest frame, Npair,reai is the number of pairs 
in an event with relative momentum fc*, and Npair,mixed 
is the number of pairs expected from random combina- 
torics and is constructed from mixed events. The region 
of interest to an HBT analysis is below k* ~ 100 MeV. 
The Npair,mixed distribution is usually normalized to the 
Npair,reai distribution at some value of k* above this re- 
gion. Our Monte Carlo model yields a correlation C{k*) 
well described by a Gaussian with an RMS width of 750 
MeV. For bubbles with ub — 8 the amplitude of the Gaus- 
sian is 8%. If the Npair,mixed distribution is normalized to 
the region near k* = 500 MeV the correlation induced by 
bubble formation would show up as a smooth 1-2% varia- 
tion of C(fc*) between and 500 MeV. In HBT analyses, 
long-range non-femptoscopic correlations of this kind are 
treated as background and ignored or parameterized. We 
therefore conclude that although HBT measurements per- 
formed at lower beam energies do not report signs of a 
strong baryon-baryon correlation as described here, they 
also do not exclude their presence. 



5. Summary and conclusion 

We have shown that bubble nucleation at a first or 
phase transition coupled with radial flow in heavy-ion col- 
lisions can lead to detectable baryon-baryon correlations. 
We've mapped out the shape of these correlations in ra- 
pidity and azimuthal angle and discussed how they will 
manifest as non-fiow or flow fluctuations by contributing 
to the difference between W2{2} and t;2{4}. We explored 
how the correlations will depend on the strength of the 
radial flow, and the temperature of the system at the time 
when the baryon rich regions decay. We reported the vari- 
ation of the azimuthal and longitudinal width with tem- 
perature, fiow, geometry, and decay time. We find that if 
bubble decay happens late in the evolution, close to the 
freeze-out when the temperature is lower and the fiow is 
larger, the correlations will be narrow in rapidity and az- 
imuth. We conclude that these correlations would be easy 
to detect. If the bubbles decay early in the evolution, the 
flow will be weaker, the temperature will be higher, and 
the correlation will be wider in both directions and there- 
fore less pronounced in the data. In this later scenario, 
there will be more hadronic rescattering which may wash 
out the signal entirely. The effect of the hadronic stage on 
these correlations is an important study which will be the 
focus of future more quantitative work. 

The observation of baryon-baryon correlations as de- 
scribed in this work will be evidence for the existence of a 
first-order phase transtion. The lack of such a signal will 



indicate that either bubble nucleation at a first order phase 
transition did not occur, or that the transition was suffi- 
ciently seperated from freeze-out as to render the correla- 
tions unobservable. We propose therefore, that studies of 
baryon-baryon correlations can be used to answer whether 
a first-order phase tranition is present in the vicinity of the 
freeze-out curve in heavy-ion collisions. Future studies of 
the effect of the hadronic rescattering stage will allow us 
to more precisely specify the region within which a first or- 
der phase transition should be detectable. If the transtion 
to quarkyonic matter is first order and close enough to 
hadronic freeze-out, then we conclude that it can be de- 
tected through baryon-baryon correlations. 
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